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Abstract 

In this paper we consider first order differential models of collective behaviors of groups 
of agents based on the mass conservation equation. Models are formulated taking the spatial 
distribution of the agents as the main unknown, expressed in terms of a probability measure 
evolving in time. We develop an existence and approximation theory of the solutions to such 
models and we show that some recently proposed models of crowd and swarm dynamics fit 
our theoretic paradigm. 

1 Introduction 

This paper deals with mathematical models of collective behaviors of groups of interacting agents, 
such as human crowds and swarms. The reference framework is that of first order differential 
models ruled by the principle of conservation of mass and supplemented by a kinematic description 
of the behavioral strategy developed by the agents. The number of individuals is assumed to be 
finite, though arbitrarily large, their state being represented by their position evolving in time in 
the Euclidean space R'', d > 1. Rather than looking at the agents singularly, we abstract their 
spatiotemporal evolution into that of a suitable probability measure, representing the law of their 
positions understood as random variables. This allows us to provide a unified theory for both 
discrete and continuous models. 

There are three main contributions of this paper, (i) We outline a basic set of modeling as- 
sumptions, which allow us to prove the existence of probability measure solutions to a broad class 
of models of collective behaviors of the kind described above, (ii) Under the very same assump- 
tions, complemented with a suitable condition on the time and space discretization, we provide a 
convergence result of an ad hoc numerical scheme, originally proposed in [3], for approximating 
the solutions to such models, (iii) We reinterpret the rendez-vous, swarm, and crowd dynamics 
models developed in [H [SJ |5J [3] in the light of our probabilistic description, showing that they 
comply with the above modeling assumptions and arc therefore in the scope of our existence and 
approximation theory. 

In more detail, the paper is organized as follows. After this Introduction, Section [01 briefly 
introduces and explains the main notations and notions used throughout the other sections. Sec- 
tion [5] proposes a probabilistic interpretation of the dynamics of systems of interacting agents, 
discussing both the indefinite mass conservation equation and the related Cauchy problems. Then 
it presents the modeling assumptions and offers an overview of the results proved in the subsequent 
sections. Section |3] deals with the existence of solutions and Section U with the convergence of 
the approximation scheme. These two sections may be skipped by readers not interested in the 
technical details of the proofs. Section [S] addresses the above-cited crowd and swarm dynamics 
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models, showing that they fit tlie theory in all cases interesting for the applications. Finally, Sec- 
tion [6] gives an example of how ODE-based discrete models can be explicitly recovered from our 
measure-theoretic framework. In addition, it shows by means of numerical tests the convergence 
of the computational scheme discussed in Section 2] to the ODE solutions of such models. 

1.1 Notations and background 

In this section we quickly review the main notations and notions that we will extensively use in 
the paper. 

Functions and function spaces. We denote by C{A; B) the space of continuous functions 
f : A B. The set B is usually omitted if it is R. Coherently, we denote by C^(R'') the space 
of real- valued infinitely differentiable functions with compact support in R"^. 

The indicator function of a set A is 1a, namely 1a{x) ~ 1 ii x ^ A and 1a{x) = if x e A^. 

Norms. We use |-| for the Euclidean norm in R'' and • for the corresponding inner product. The 
open Euclidean ball with center a; S R'' and radius i? > is denoted by Bii{x). If / : R'' ^ R'' is 
bounded, its oD-norm is ||/|loo := sup^jgRd |/(a;)|. 

Measures. Let B{R'^) be the Borel cr-algebra on R*^. If /i is a measure on the measurable space 
(R'^, B{R.'^)) and / : R'' ^ R'^ is Borel, the integral of / w.r.t. /i over a measurable set A is denoted 
by f dfjL, or by f{x) d^{x) when it is necessary to emphasize the variable of integration. The 
Lebesgue measure in R'' is denoted C^. However, for integrals with respect to we will prefer 
the usual symbol dx to dC'^{x). If / is Borel, we denote by the push forward of fi through 
/. Specifically, is the measure defined by the relation 

J 77d(./#/i) = Jvofdfi (1) 

for all bounded (or /#/i-integrable) and Borel t] : M.'^ ^ M.. Taking rj = 1a, A E B{M.'^), gives in 
particular = ^lir\A)). 

Probability spaces. We denote by 7'i(R'^) the space of probability measures on (R'^, S(R'^)) 
whose first moment is finite, i.e., J^a\x\d^{x) < +oo. The space 7^2(R''), that we will also 
occasionally mention, is defined analogously using the second moment. Given two probability 
measures v € 7'i(R'^), their Wasserstein distance is defined to be 

Wi{iJ.,v)^ sup tpd{i^-fj.), 

where Lip2(R'') is the space of Lipschitz continuous functions (yS : R'' R with Lipschitz constant 
LipCiys) < 1. It can be shown that Wi is a metric on 'Pi{M.'^) and that (7'i(R''), Wi) is complete 
(see e.g., [U Proposition 7.1.5]). 

Finally, to deal with curves in 7'i(R'*) parameterized by time, [0, T] 3 t i-^ fit £ Vi{M.'^), we 
introduce the space C([0, T]; 7'i(R'')), which is complete with the metric 

dist(/x,, u,) sup Wi{fit, vt). 

*e[o,T] 

2 Problem statement and main results 

In this section we present our approach to the modeling of systems of interacting agents by means 
of probability measures, and we give an overview of our results. 
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2.1 Probabilistic description of systems of interacting agents 

Crowds and swarms can be thought of, in abstract, as systems of N interacting agents in the 
physical space M.'^. The evolution of such systems in time is described by tracing the agent positions 
at successive instants. Assume that the position of the i-th agent at time i is a random variable 
XI from a fixed (i.e., time-independent) abstract probability space {il, P) to the measurable 
space (K"^, B{W^)). The probability P is naturally transported by XI onto the new probability 
^^ ■- Xl^P on {W^, B{W^)), called the law of XI The fact that nt does not depend on the 
agent label i means that agents are indistinguishable: the probability of finding a certain agent 
somewhere in at time t is the same for all agents, namely, given A E B{M.'^), P{{Xl G A}) = 
fitiA) for aWi^l, N. 

Let us now fix a measurable set A C and count the average number of agents contained in 
A at time t. This amounts to introducing the new random variable Yt^A : — >■ N defined as 

N 

Yt,A = '^{xieA} 

i=l 

and taking its expectation, that we can compute as follows: 

N 

E[Yt,A] - e A}) = N^^M)■ 

i=l 

Notice that E[yt_,], thought of as a map on i3(R'^), is a finite positive measure, say mt, over the 
measurable space (M'', B{M.'^)), such that m,t{W^) = N all t. It is straightforward to identify mt 
with the mass of the system at time t. From the above calculation we see that m* is proportional 
to the probability measure of the distribution of the agents: mt{A) = N^t{A), where mt{A) is 
the mass of the set A E B{W^) at time t. In addition, mt(M'^) = N for all t, i.e., the total mass of 
the system is constant in time. 

The latter observation suggests that we can assume the principle of conservation of the mass, 
stating that the mass of any measurable set A may change in time only because of infiow or outfiow 
of mass from the boundary dA. In other words, the mass is neither created nor destroyed but 
only moved across the domain. This is expressed by postulating the continuity equation (or mass 
conservation equation) for the evolution of the measure m^: 

^ + V • {mtvt) = 0, (2) 

where Vt{x) is the velocity at time t in the point x e M''. In systems of interacting agents 
the velocity is likely to be affected by the distribution of the agents themselves. Due to the 
proportionality between mt and fit , this implies that Vt rnay ultimately depend on the probability 
/it. We assume in particular that vt depends on t only through fit itself, i.e., that the system is 
autonomous. Finally, we write Vt = v[fj.t] to emphasize such a structure of the velocity and notice 
that Eq. ^ can be converted into an evolution equation for the probability fit '■ 

^ + V.ifitv[fit]) = (3) 

for a; € M'' and t G (0, T], where T > is the final time. If fit solves Eq. Q, then mt formally 
solves Eq. ^ with vt = v[mt/N]. 

2.2 Cauchy problems 

By supplementing Eq. with an initial condition /2, the following Cauchy problem is obtained: 

■{fitv[fit])^0 inM'^x(0, T] 
fio — ft in R'', 
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which models the spatiotemporal evolution of the agent distribution starting from the initial 
configuration described by fl. Derivatives in problem (|4|) are meant in the sense of distributions, 
which leads us to consider the following notion of weak solution: 

Definition 1 (Weak solutions). Given p, G 7'i(M'*), we say that (x, e C([0, T]\ ViiW^)) is a weak 
solution to problem ^ if 

t 

J (j)dfit= J (t)dfl + J J v[fir]-V(t)dfirdT, V0e C,°°(R''), Vt e [0, T]. (5) 

R<i W Rd 

As far as the probabilistic interpretation is concerned, we take the initial condition into account 
by understanding fit in as the law of the random variable XI :~ K[XI\Xq], the expectation of 
XI conditioned to its initial value Xq. In practice, fit is reinterpreted as the distribution of the 
position of the i-th agent subject to the distribution of the corresponding initial position. From 
the theory of conditional expectation, we known that X^ is a function of Xq, i.e., there exists a 
Borel mapping 74 : M'' ^ such that X^ = ^t{Xl). This implies, in particular, 

Mt = 7t#A- (6) 

In order for Eq. (|6]) to result in a representation formula for the solutions of problem (|4]), a 
more precise characterization of the function 74 is needed. Formally, we plug Eq. (jB]) into Eq. ([5]) 
and compute 

t 

[<l}{lt{x)) - (j){x)]dfi{x) =11 v[-1r#lA{lT{x)) ■ S/(j){-frix))dfl{x)dT 



Rd 



for an arbitrarily fixed test function Next we notice that the integrand at the right-hand side 
can be read as the derivative w.r.t. t of the function 0(7T-(a;)), provided we identify ■§tJt{x) with 
z;[7t#/2](7((a;)). Under this assumption we get 

[(l^iltix)) ~ 4>{x)] dfl{x) ^ J J -^Hlrix)) dfl{x) dr 

Rd Rd 

and further, interchanging the order of integration at the right-hand side, 

[(j){jt{x)) - (j){jo{x))]dp{x). 

With the additional condition 70(2:) = x (i.e., 70 is the identity function in Mf^), this shows that 
fit represented by Eq. ^ is formally a weak solution to the Cauchy problem 
To sum up, "ft has been characterized as a function such that 

''^^ = t'b#A](7t(^)), te{0,T] 

Jo{x) = X 

for every x GMf^. In transport theory, such a function is called a flow map. The physical interpre- 
tation is that 7t(A) is the configuration assumed by the set A at time t > when transported by 
the velocity field v. Alternatively, for x S M'' the mapping t i-^ jt{x) is the trajectory of system 
([7]) issuing from x. 

The method of representing solutions to the Cauchy problem (|4]) via flow maps is called method 
of the characteristics. To develop our theory we will mostly prefer a different approach, more 
suited to treat, by common ideas, existence and approximation of solutions to the models we are 
interested in. The reader interested in the method of the characteristics is referred to [5] and 
references therein for further details. 
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2.3 Basic assumptions and results overview 

Models based on Eq. ([3]) require to specify the velocity, namely its dependence on the probability 
fit and on the space variable x. Rather than considering a specific model, we outline here a small 
set of assumptions, which can be possibly regarded as modeling guidelines, whence the whole 
theory will follow. 

Assumption 1 (Properties of v). We assume that the velocity field v = v[fi]{x) satisfies the 
following properties. 

(i) Uniform boundedness: there exists V > such that 



(ii) Lipschitz continuity: there exists a constant Lip(w) > 0, independent of both the space 
variable and the probability, such that 

\v[iy]{y) - v[fi]{x)\ < Upiv){\y - x\ + Wi{fi, ly)), Vx, y G V^, e ViiR"). 

(iii) Linearity w.r.t. the measure for convex combinations: 

v[afi + (1 - a)iy] = av[fj] + (1 - a)v[iy], \/ fi, v e 7'i(R''), Va G [0, 1]. 

From this basic set of hypotheses we will be able to prove, in Section [31 existence of solutions 
to the Cauchy problem ([4]) in the appropriate weak sense of Definition [T] More precisely, to this 
end we need the further technical assumption that the initial condition have finite first and second 
order moments: 

Assumption 2 (Initial condition). We assume that fi e Vl[m.'^)c^V2{^'^)■ 

Then, our main result in Section [3] reads: 

Theorem (cf. Thcorem[5]). Under Assumptions{^\^ there exists a weak solution to problem ([¥]). 

Notice that Assumption [5] is readily satisfied if, for instance, fl has compact support. Indeed, 
in such a case supp /I is bounded, i.e., there exists a ball Bj^ {Q) of sufficiently large radius R > 
such that, for every p > 0, 



The compactness of the support of /2 makes perfectly sense from the modeling point of view, in 
fact a crowd or a swarm spread on the whole space would sound quite unrealistic. Assumption [2] 
is therefore not restrictive for our purposes. 

In Section |4] we turn our attention to the approximation of solutions to problem (|4]) . We 
introduce a sequence of grids in R'' x [0, T] with mesh parameters h^ in space and At^ in time. 
The index k relates to the grid refinement, in such a way that hk, Atk — > when fc — !■ oo. 
Specifically, we consider the numerical scheme proposed in [S], which at each time step seeks 
an approximation of /Zt via a probability measure absolutely continuous w.r.t. to Lebesgue and 
piecewise constant in space. By introducing a linear-in-time interpolation of such approximate 
solutions and passing to the limit fc — > oo, we obtain the following convergence result: 

Theorem (cf. Theorem [T^ . Under Assumptions[ll[^ suppose that hk = o{Atk) for fc oo. // 
the sequence of approximate solutions converges to some n, S C([0, T]; 7'i(R'*)) when the grid is 
refined then ji, is a weak solution to problem 



\vWx)\ <V, Vx e R'', V^t e 7'i(R''). 
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Note that convergence when the grid is refined is an assumption of this theorem. In this respect, 
this result resembles the Lax-Wendroff 's Theorem for the numerical approximation of hyperbolic 
conservation laws (see e.g., [?])• However, we anticipate that if there exists a bounded subset of 
M.'^, that at each time step and for all level of refinement of the grid contains the supports of all 
the approximate solutions, then the sequence docs converge in C([0, T]; 7^i(K'')) to some limit, 
which is then a weak solution to problem (jH) (cf. Corollarv [T5|) . 

In Section [5] we apply the above theory to the models of swarm and crowd dynamics presented 
in [SJ IH] . As shown in [5] , these models can be derived from the common framework provided by 
Eq. ^ , with a velocity field of the form 

v[fit]ix) = Vdix) + N J fi\y - x\)r{y - x)xuAy)dfJ.tiy), 

where the integral expresses the interactions among the agents. Assuming some minimal regularity 
of the functions v^, /, r, xu^: we prove that this velocity complies with Assumption [1] in case of 
both isotropic and anisotropic interactions, cf. Sections l5.1l and [5^ respectively. As a consequence, 
for the above-mentioned models we deduce existence of probability measure solutions, that can be 
duly approximated via the numerical scheme discussed in Sectional 

Our results are exemplified in Sectional where we show that, given a purely atomic initial 
measure /2, i.e., 

1 ^ 

1=1 

a solution to problem Q can be found by solving a system of ODEs whose unknowns are the 
trajectories of the agents. In addition, using a numerical solution of these ODEs as a benchmark, 
we are able to visualize the convergence of the numerical scheme presented in Section |4l 



3 Existence of solutions 

This section is devoted to give a constructive proof of the existence of solutions to the Cauchy 
problem Q. Under Assumptions [TIIH the solution is constructed as the limit of a suitable sequence 
of curves in ViQS.'^) parameterized by time t. 

Let {Atk)k>o be a sequence of time steps such that Atk when k — > oo. We consider the 
measures (/^^)n>o generated recursively as 

fA^Ji+i =7^#mJ;, n = 0, 1, ...,iVfc-l, 

where Nk € N is such that NkAtk = T and 7,^ is the (one-step) flow map 

j'^{x)^x + v[f,'^]{x)At,. (9) 

It can be shown (see e.g., [9]) that ([8]) is the explicit time discretization of (HJ at the time instants 
t!^ = nAtk- In particular, the mapping 7^ results from the explicit time discretization of problem 
([7]). Since /I is a probability measure, by induction it is immediate to check that so are all of the 



By linear interpolation in time, we define the following curves: 



n=0 



At, At, ^"+1 



Obviously, Af^ is a probability measure for each t £ [0, T] and each A; > 0. 
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We will use the curves to construct, in the limit k ^ oo, a weak solution to problem (j?]). 
The proof is divided in two parts, each of which proceeds through a series of technical intermediate 
steps developed in the next two sections. First, wc show that the measures Mj^ converge to a limit; 
later, we show that such limit satisfies Eq. ([5]). 

3.1 Convergence of the measures 

In this section we prove that, when fc -> oo and up to subsequences, converges to a limit /i, 
in C([0, T];Vi{R'')). 

We begin by establishing the necessary regularity properties of the iterates and the curves 

Lemma 1. We have fi^ £ 'Pi(R'^) fM^2{^'^) for all n = I, . . . , Nk and allk>0. In addition: 

sup sup / |a;|^ d^^lx) < +oo, p = 1, 2, 

fe>0 l<ri<JVfc J 

i.e., first and second moments of the fi^ 's are uniformly bounded. 
Proof. (i) Let us begin by considering the case p = 1. Note that 



< J \x\ dn'^{x) + Atk J \v[fit]{x)\ diil{x) 

< j \x\ dfx'^ix) + VAtk, 



then 



\x\ / \x\ df4ix)<VAt,, 



and summing telescopically over n we get 

\x\ dfi'^ix) < I \x\ dp{x) + VT. (10) 



Since ft £'Pi(R''), this implies /ijj e 7'i(R'^) for all n = 1, . . . , Nk and aU A: > 0. But the 
right-hand side of (jTP)) is independent of both n and fc, hence this also provides the uniform 
bound for p = 1. 

(ii) We argue analogously for the case p = 2. Wc note that 

dfj.'^+iix) = / \j^{x)\^ d^l'^{x) = / \x + v[^l'^]{x)Atk\^ d^l'^{x) 



and that 

\x + v[^i'^]ix)Atk\^ = \x\' + 2AtkX ■ v[fi'^]ix) + At'i 
< \x\^ + 2VAtk\x\ + V^Atl 

therefore 

d^^'^^,{x) < J \x\' dti'^{x) + 2VAtk J \x\ dfi'^ix) + Atl 
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Moreover, using the bound ([TU)) we deduce 

J \x\^ - J \xf dMf;(^) < 2FAtfc J \x\ dfji{x) + 2V''TMk + 

and, summing telescopically over n, 

\x\^ d^il{x) < I \x\^ dfL{x) + 2VT I \x\ dfi{x)+3V^T^ 



whence the claims of the lemma follow also for p = 2. □ 
Lemma 2. For all k > we have 

Wi{M^,M^)<V\t-s\, Vs, i€[0, T], (11) 
i.e., the curves Mj^ E C!{[0, T]; 7'i(M'^)) are Lipschitz continuous uniformly in k. Moreover, 



sup sup / |a;|P dM^ix) < +oo, P = 1, 2. 



ne[o,T] 

Proof. (i) Wc claim G ViiR'') for all t, k. To show this, wc fix fc > and t G [0, T] and 
observe that there exists < n < Nk such that t € [t^, i^+i]. Hence, using Eq. (fTOj). we 
obtain 

J \x\ dM^{x) ^[l- ^) j \x\ d„'^{x) + ^ J \x\ < J \x\ dUix) + VT. 

From the arbitrariness of t G [0, T], fc > our claim follows, along with the uniform bound- 
edness of the first moment of Mj^ in both t and k. 

(ii) We prove now the estimate ([TT|) . Let s, t G [0, T] and assume, without loss of generality, 
that s < t. Then there exist two integers m and n, such that < m < n < Nk, with the 
property that < s < t^_^i and t^^ <t < tjj^i. Therefore we can write 

rk _ f -I * ^ tm \ ,,k , ^ ^ ^rn ,,k l\ fk _ f -i ^ ^ ,,k , ^ ^ ,k 



and further, owing to the triangle inequality, 

n-l 

VFl(A^^Mf)<VFl(M^M™+l)+ E W^i(a*J,A*'+i)+W^i(m', M')- (12) 
Notice that 

W, {Ml fit^,) ^ sup / ^ d{fit+, - Ml 

V3eLipi(R'') J 



R"* 



At_ 

and analogously Wi{fi'^, M^) = {t - t'^)Wi{^';^, /i^+i)/Atfc, thus, according to Eq. ([T^ . 
estimating H^i(Mg'^, M/^) amounts to estimating Wiinl Hi+i) for arbitrary < i < A^fc — 1. 
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For any ip £ Lipj(M''), using = lii^fJ-i yields 

fix) d{fi'^+i - ti'i){x) = / ifilH^)) - fix)) d^i^ix) 



< 



j |7f(x) - x\ dfi'^ix) = AtkJ \v[tJ,'i]{x)\ dtx,{x) < VAtk, 



consequently Wi{^^, Mi+i) ^ VAtk all i. From Eq. wc deduce 

t¥l(A^^ M^) < V[{C+, -s) + {n-m- l)Ah + (< - t)] 
- F(t - ,s + - 4 - Atk) = F(i - ,s), 

which proves our claim, 
(iii) Finally, we claim Mj^ £ V2{^'^) for all t, k. Arguing like in (P we have 

t — t'^,\ f . 2 ^,,k(\ , ^ ~ [ I |2 



< J \xf dfL{x) + 2VT J \x\ dfl{x) + 3V^T^ , 

where Lemma [1] has been used. This proves the claim and also the uniform boundedness of 
the second moment of in both t and k. □ 

We are now ready to prove the main convergence result of this section. 

Proposition 3. There exists /i, £ C([0, T]; Vi{M.'^)) and there exists a subsequence 7">0 
such that 

lim sup Wi{Mt', ^it)=0. 
J^°°ie[o,T] 

Proof. It suffices to prove that {A/^}fc>o is a relatively compact subset of C([0, T]; 'Pi{M.'^)). Owing 
to Ascoli-Arzcla's Theorem, this happens if {il/^}fe>o is equicontinuous and {M^}k>a is relatively 
compact in ViiR'^) for aU t £ [0, T]. 

(i) Equicontinuity follows from the estimate ([TT|) . Indeed, let e > then for S = e/{2V), which 
does not depend on k, we have |i - s| < (5 Wi{M^, Mj^) < e/2 < e. 

(ii) According to [H Proposition 7.1.5], the relative compactness of {Af/^}fe>o in ViiW^) is equiv- 
alent to the fact that {M^}k>{) be tight and have uniformly integrable first moments. 

(ii-a) Using [1] Remark 5.1.5], a sufficient condition for tightness is that there exists a function 
i^p -.W^ ^ [0, +00], whose sublevel sets {x e M'' : Lp{x) < c} are compact in M'^, such 
that 

sup / (fidM^ < +00. 

fe>0 J 

Taking ip{x) = and invoking Lemma [5] we see that this condition is fulfilled, hence 
{M/^}fc>o is tight. 

(ii-b) Using [1] Eq. (5.1.20)], a sufficient condition for the uniform intcgrability of the first 
moments of {M^}k>o is that there exists p > 1 such that 



sup / dM^{x) < +00. 

fc>0 J 



From Lemma m we know that this actually holds for p = 2. 

Since {M^}k>o is relatively compact in C([0, T]; Vi{M.'^)), up to subsequences we obtain that the 
sequence (M^)fe>o converges in C([0, T]; 7'i(R'')) and we are done. □ 
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3.2 The limit /x, solves problem (jl]) 

In the previous section we have constructed a map /i, as the hmit of the sequence {M^)k>o- In 
this section we prove that such a /x, is a weak solution to problem To this end, we first derive 
an equation solved by the Af/^'s, then we pass to the limit fc — > oo thanks to Proposition [31 

In the following, (p e C'^iW'-) is a fixed test function. Using the definition of M^ , let us 
compute 



n=0 
Nk-1 

= E 

n=0 



1 - 



t-t 



Atk 



(13) 



A Taylor expansion of </) o with Lagrange's reminder gives 

<^(7^(x)) =0(x + w[Mji](a:)Atfc) 



0(x) + V0(x) • v[^li]{x)Atk + -{D^{x)v[^il]{x)) ■ v[^ii]{x)At 



where is the Hessian of (j) and a; is a point of the segment connecting x and x + v[iJ^{x)Atk- 
Hence the previous computation specializes as 



^dy^l + it-tl) / i;[M^].V0dM, 



-Atk{t-t':,) / (Z?2^(x)7;[AiJ;]) -i^/i'lrfM' ) 



We claim now that the mapping 1 1— > Jj^^ (/> dM/^ is Lipschitz continuous, hence a.e. differentiable 
by Rademacher's Theorem. To see this, observe that x i-> 0(x')/ Lip((^) is Lipschitz continuous 
with Lipschitz constant at most 1, so that 



(pdMl 



= Lip(0) 



Lip((/.) 



d{M^ - M^) 



< Lip(0)I^l(M,^ M,') < Up{cb)V \t-s\. 
Thus, using Eq. p3p . we compute the derivative 

Nk-l 



dt 



■,dlvL 



R<i 



n=0 



Let us consider now: 



Atk 



v[M,^]■Vcpd^^'^^Al[,.^,.^^^it) 
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and invoke Assumption \T[ (pii)) to get 

Nk-l 



J2 W^M^lK,t:;+jW 



Using this in ([14)) gives 



t - 



t - 1 



At 



k 



At 



(15) 



Formally we can regard this expression as an equation satisfied by M/^. The remaining of this 
section is devoted to relate this equation to Eq. ([5]). 

Set H := ||Z?^0||oc = Lip(V0) and notice that, owing to Assumption [I]- , the function 
X I— >■ u[^j!;](a;) • V(/)(a;) is Lipschitz continuous with 

Lip(iJ • V0) <HV + Lip(t;)||V(/)|loc^ L. (16) 

Consequently: 



(j)dM^ - / v[M^] ■ \/(f>dMt 



n=0 { 



t ~ f'' 



2 



Atk 

2 



Atk 



<cY.\Atu + ^-^W,{f,lf,t+i) 



n=0 
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where C := m'Ax{^HV'^, L, 2\\S7(p\\oc}- But Aifj+i) < VAtk, therefore 



< 



C At, + V{t-t'^) + V^-—^ 1^, ,.^^] 

n=0 ^ k ) 



it)- 



Now, for t e [t^, it rcsuhs t - t'^ < At,. In addition, X;!^=o ^ ^[tt tJi+J (*) = 1[o,t](0: hence 

finahy 

<C{l + 2V)Atkl[o,T]{t)- 
Integrating Eq. between and t < T and using the last inequahty yields 



\ M<i 



dT 



< C(l + 2V)TAtk 

so that, by further manipulating the left-hand side and taking the limit for fc — > cx), we obtain 



lim 



Rd 



0. 



(17) 



To infer from (flT)) that /^t solves (O, we need the following convergence result. 
Lemma 4. When k oo we have, up to subsequences, 

t t 
J (|)dM^ J (pd^it and J J v[M';] ■ V(f>dM^ dr ^ J J v[fir] ■ (l> d^ir dr 











for allte [0, T]. 

Proof. (i) For the first limit we write 



J (l)dfit - I (t)dM^ 



^d{pt-M^) 



Lip(0) 



Lip(0) 
<Lip((/.)iyi(Af^A^O, 

and, up to passing to a suitable subsequence of (M^)fc>o, we conclude by applying Proposi- 
tion [31 

(ii) For the second limit we preliminarily observe that | v[M^] ■ V^dM^'j < V^||V0||oo, thus, 
by dominated convergence, 



Jim j j v[M^] ■ \7(l)dM^dT = J i \im J v[M^] ■ \/(f>dM'^ J dr. 



(18) 
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Now 



< (||V0||ooLip(«) +i)W^i(Af^ ^lr), 

L being the constant defined in Eq. Proposition [3] implies then J^^ v[M^] ■ S/(f>dM^ — > 

Jjgj v[iIt] ■ V4)d^ir for all r € [0, T], and the thesis follows from Eq. ([TS]) . □ 

Combining Eq. (|17p and Lemma 21 and thanks to the arbitrariness of 0, we obtain that /i, 
solves problem Q in the sense of Definition [TJ In conclusion, we have proved: 

Theorem 5 (Existence). Let Assumptions hold. Then there exists a weak solution /it, G 
C([0, T]; 7^1 (M'')) io i/ie Cauchy problem g]). 



4 Approximation of the solutions 

This section is devoted to a convergence analysis of the numerical scheme proposed in [3] for 
the approximation of the solutions to problem We begin by sketching the main ideas which 
underlie the construction of the scheme, referring the reader to the above-cited paper for a detailed 
derivation. The scheme is obtained from a twofold approximation of problem (0]) , in time and in 
space. To this goal, we introduce a discretization of the time interval [0, T] by means of discrete 
instants t^ = nAtk, where the index n ranges from to a value Nk such that Nk/^tk = T and the 
time step At^ > tends to when fc — > cxj. By this discretization, we obtain the discrete-time 
dynamical system Q. Then, we introduce a space discretization in the following way. We define a 



pairwise disjoint partition of R'^ made of measurable elements E B{1 



in 



Id) G 



fc > 0, such that 



u 



for all A: > 0. For the sake of simplicity, we assume that the E^'s are hypercubes of edge length 
hk > 0, with /ife — > when k — > oo. Specifically, 



Ef 



d 

X 

1=1 



Note that the index k identifies the level of refinement of the numerical grid. Using this space 
discretization, we approximate both the initial condition p, and the measures fJ^ by means of 
piecewise constant measures Aj; ^ , n = 0, . . . , N^, k > 0. More precisely, dX'^ = pj^ dx with 



The space discretization makes it necessary to approximate also the flow map -f^^ deflned in Eq. 
This is accomplished by the mapping 



^^{x)=x + v^^{x)Atk 



it 



ix) = J2v[X^„]ix'l)lE.{x), 



x^ being a point of the grid ceh E^ (e.g., its center x'^ = ihk). In practice, the velocity z;[//J;](a;) 
is approximated by the piecewise constant field ^^(a;) taking in E^ the value that v, computed 
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w.r.t. the measure A^, takes in . Notice that < V for aU x G R'', aU ?i = 0, . . . , Nk, and 

all fc > because of Assumption [TJ-dll)- 

By imposing Ajj_|_j(i?f) — (7,^#Aj;)(i?f ) for each i g Z'^, we deduce the following explicit-in- 
time scheme relating recursively the coefficients {p"} at two successive time steps: 

p"^' = F E p^^(^' ^ (inrHEt)), (19) 

k jg^d 

where hf is C'^{Ef). To start up the scheme one has to provide the coefficients p", that we obtain 
as 

p^ = ^/2(i?f), tez^. (20) 

We are now ready to start the analysis of the proposed scheme. To this end, we begin with a 
simple property of the approximation measures. 

Lemma 6. For every n = 0, . . . , N^, k > 0, the Af^ 's are probability measures. 

Proof. (i) The claim is certainly true for Aq, since Eq. ((20)) shows that p° > all i and further- 
more 

iez<i tel.'' 

where we have used the cr-additivity of /2. 

(ii) If we assume now that A^j is a probability measure for a certain n, using Eq. ([T^ we get 
p^+i > alH and moreover 



Af,+i(R'^) = hi p^' = E E ^"^'(4 n iinr\E^)) 

= E ^" E ^'(4 n (7^)"'(i?- )) = E ^"^'(4 n (7^)"'(K'^)) 



= /.^E^" = ^"(^') = i' 

where we have used the fact that (7,^)~"'^(£'f ) are pairwise disjoint and that set operations, 
e.g., union, commute with the inverse image of a function. By induction on n, and by the 
arbitrariness of fc > 0, the claim follows. □ 

The next result shows that the link between two successive measures Afj, Aj^^^ is a push forward 
as defined by Eq. ([1]) , provided we restrict test functions to simple functions adapted to the spatial 
grid {E^}ifzzd, i.e., pieccwise constant functions s -.W^^M. of the form 

s{x) = E] Q!il£;fc(a;) {cxi G M). 

iei,'' 

Lemma 7. Let s : — > M 6e simple over the grid {E^ }j^^^d. Then 
Proof. We have 
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Notice that s(7^(x)) = J^iez-^ '^i'^{^'')-^(.E''){^)j i-^-i s o is piecewise constant over the sets 
{j^)~^{E^), which form a pairwise disjoint partition of R'^. Therefore, for any given A £ B{Mf'-), 

For A = this enables us to continue the previous computation as 



and to obtain the thesis. 



□ 



An immediate consequence of this lemma is the following result, that will be fundamental for 
the sequel. 

Lemma 8. For all Lipschitz continuous have 



<2Up{ip)Vdhk, Vti = 0, . . . , iVfe. 



Proof. We consider the simple function 

six) = Vi^'^)'^E^i^ 

and, using Lemma [71 we compute: 



^(a;)dA^+i(x) = / ((^(x) - s(a;)) dA^+i + / s{x) dX'^^,{x) 



Mx) - s{x)) dXi+,{x) + / [s(7^(x)) - <p(7^(a;))] d\l{x) 



MA 



+ ^(7^(x))dA:;(a;) 



Y\ {v{x)-^{x^))d\i^^{x)+ /(^(x,^)-^(x))da';#A,';)(^) 



+ ^(7:(x))dAf,(a;). 



Notice that 7„#A„ is a probability measure. The Lipschitz continuity of ip entails 



ip{x)d\l+^{x)- / ip{^l{x))d\i{x) 



x~xf\ dX'^,^{x)+ / \xf-x\ d(7*#A^)(a;) 
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whence, since Ix — I < diamSf = Vdhk for all x G E^, 



and, due to the arbitrariness of tp and n, the thesis follows. 



□ 



To address the convergence of the measures A^, it is convenient to introduce the following 
linear-in-time interpolation: 



A 



Nk-l 

E 

n=0 



1 - 



f _ fk 

''71 

Atfc 



'' ''n \k 



(21) 



with a probability measure for all t e [0, T] and all k > 0. This definition is reminiscent of the 
definition of M/= (cf. Section [3l). 

In order to prove a convergence result about Aj', from now on we make the following assumption: 

Assumption 3. The mesh parameters hk, Atk satisfy 

hk = o{Atk) for fc — > oo. 

As a consequence, there exists a sequence (/3fc)fc>o C K, with /3k > all k and limfc_>oo l^k ~ 0, 
such that hk = /3feAtfe. By convergence, it further results supj.>Q h^, supj.>Q At^, supj,>Q < -|-oo. 

4.1 Regularity of 

This section is devoted to establish some regularity properties of the curves A^. 
Lemma 9. We have A^j G 7'i(R'') for all n = 0, . . . , Nk and all k > 0. In particular, 

sup sup / |x| dX'^{x) < +00, 

fc>0 a<n<Nk J 

w 

i.e., first order moments of the Ajj 's are uniformly bounded. 

Proof. (i) We begin by proving an inequality for n = 0. Recalling Eq. (|20p . we have 
J \x\ dX^,ix) =Y.P'J 1^1 '^•^ = E p / 1^1 dxfl{Ef) = J six)dflix), 

where s is the simple function s = X^ieZ"* with a; = J^k \x\ dx. Thus 



|a;| dX^ix) - J \x\ d^i{x) 



< J \s{x) ~ \x\\ dii{x) 



Considering that \ai - |a;|| = /gfedyl - la^Dc'yl < 77- Je^ N - 2/1 ^2/i we further deduce 

-Ep/ J \x - y\ dydfi{x) < Vdhk, 



E'' B * 
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where we have used that, in the double integral, both x and y are points of the same grid 
cell E'f , hence |x — y| < dianiE'f = Vdhk- It follows 

|a;| dX^ix) < / |x| dfl{x) + Vdhk, 



whence the thesis for n ~ 0. 

(ii) We obtain the general case < n < by induction from this and Lemma \8\ Choosing 
f{x) = \x\ in the latter yields 



\x\ < j |7^(x)| d\^^ix)+2Vdhk 

J \x\ d\i{x)+^tk j dX'^ix)+2Vdhk 
J \x\ d\l{x) + VMk + 2Vdhk 



< 

< 

whence, summing tclcscopically and using ~ l^kAtk, 



\x\ dXl{x) < j \x\ dX'^ix) + n{VAtk + 2Vdhk) 

< / dpix) + {V + 2Vdl3k)tn + Vdhk. 



The thesis now follows taking the supremum of both sides in n and k while considering that 

ti<T. □ 

Lemma 10. For all k >0, 

Wi{A';,A'^)<{V + 2Vd^)\t-s\, (22) 

where (3 := supk>o l^k> curves e C!{[Q, T]; Vi{W^)) are Lipschitz continuous uniformly 

in k. Moreover, 



sup sup / \x\ dA^{x) < +O0. 
k>0 te[o,T] J 



Proof. (i) Fix fc > and t G [0, T]. There exists < n < TV^ such that t e [t„, i„+i], hence 
\x\ dA^ix) ^fl- I \x\ d\l{x) + ^ / \x\ d\i^,{x). 



Owing to Lemma [HI we can find a uniform upper bound on the first moments of the A^'s: 
there exists a constant C > 0, independent of n and fc, such that 

which says that G ViiW'-) for all fc > and all t e [0, T] with uniformly bounded first 
moment. 
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(ii) We argue the Lipschitz continuity of the mapping i i— > as in Lemma [5] In particular, for 
s, t € [0, T] with s < t we obtain 

fk _ f _ fk 

W,{A':,A'^)<^I^^^W,{X'^,X'^^,)+ VFi(A^^ A^\0 + ^W^i(A,';, Af^^O, (23) 

J— m+l 

where < m < n < Nk are such that s £ [t'^, i^+i] ^-iid t E In order to estimate 

VFi(Af , Xi^i) for a generic < i < A^fe we fix <^ € Lipj(K'') and use Lemma |SJ 

ipdiXl, - ) < / i^itix)) - H^)] dXlix) + 2V~dhk 



j |7^(x)-x| dX'^{x)+2Vdhk = Atk J \v^{x)\ dX'^{x) + 2Vdhk 



< 

< VAtk + 2Vdhk ^{V + 2VdPk)Atk. 

Consequently W^i(Af , X^j^^) < CAtk all i, where C -.^ V + 2\/dsupj.>o /3fc- Computing as in 
LemmalU in view of Eq. (^5)) this yields Wi{A'^, A^) < C{t — s) and we have the thesis. □ 

4.2 Limit equation for 

The next step is to find, similarly to what we did in Section 13.21 an equation satisfied by A^' in 
which to pass to the limit fc — > oo. Fix a test function € C^(R'') and notice that the mapping 
t H- > Jjjd (j) dA(' is Lipschitz continuous because so is the mapping t i— > A^ in view of Lemma 1101 
Then, owing to Rademacher's Theorem, it is a.e. differentiable. Using expression (|2ip . we find 
that its derivative is 

d f 1 ^''^^ f 

Jt J ^ATk^ J ^'^^^'+' ~ W- (24) 

Let us introduce now the function : M'^ M'', 

9nix)^x + v[X'^]{x)Atk, 
i.e., the flow map 7^ (cf. Eq. ^) computed w.r.t. the measure Aj; instead of /x^. Then 

whence, expanding (f>o in the second term at the right-hand side with Lagrange's reminder (cf. 
the analogous calculation performed in Section [3T2] for the function <j>oj^), 

R<* J 

E |/«[An]-V</>dA:; + iAi,.|(i^2^(s>[A:;])-«[A:;]dA:;U 



At 

Nk-1 



n—0 
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where x is a point of the segment connecting x and x + w[A^] At^ 



On the other hand, computing as in Section [3.21 we find 

Nk-l 



t - r 
Si 

Atk 



hence finahy 



dt 



^ I -At, j{D^mv[xi])-v[xi]d\i 



n=0 

R 

Atk 
Atk 



v[X>^]-Vcbd{Xt+,^Xt) 

{v[x^,+,]^v[xt])-y^dxt 



This is formaUy an equation satisfied by A,. A derivation analogous to that performed in Sec- 
tion 13.21 gives 



dt 



0dAf - / i)[Af] • V(/)dA 



R-i 



< 



At 



n=0 



0rfA^;+l - / 0og,^dA^ 



+ {Atfc + Wi(A^ Af,+i)} Ik.*:;^,] W> 



n=0 



where C > is a constant depending only on V, Lip(w), Lip(V0) = HD^^jloo- In addition, from 
Eq. (EH) we know Wi{X'^, A^+j) < (V + 2Vd$)Atk, thus 



Ni.-l 



< 



At 



1 

7:T. 



n=0 



^dX';,^, - I (^og^^dXi 



C being a new constant including the previous one and V + 2v^/3. 

To estimate the remaining term at the right-hand side, we need to adapt Lemma [5] to g^. 
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Lemma 11. For all Lipschitz continuous ip 



have 



V dXn+i - V° 9n dK 



< C Lip {ip)\/dhk 



where C := 2 + Lip(w) supj,>Q Atk- 
Proof. It suffices to observe that 



AtkUp{ip)Up{v)J2 / - x\ dX'^.ix) 



< VdhkAtk Lip(iy9) \Ap{v). 
The thesis follows from this inequality, combining the triangle inequality and Lemma [S] 
With this result, we can further manipulate the previous inequality and obtain 



□ 



< C(/3fc + Atk) 



which, integrating both sides in time from to t < T, implies 



/ ( ^ / 'I'dA'; - J v[A';] • V0dAM dr < J J (j)dA^- J v[A^] ■ S/cjjdA'; 



dT 



< CT{l3k + Atk 



and finally, taking the limit fc — >■ cxj. 



lim 



0dA« - / 0dAS - 



J J v[A^]-V(l)dA'; dr 



0. 



(25) 



Thanks to Eq. (j25p . we are in a position to prove our convergence result of the numerical 
scheme. 

Theorem 12. Let Assumptions[l}^hold and assume that the sequence (A^)fc>o converges to some 
/i, in C([0, T]; Vi{M.'^)) when k oo. Then /i, is a weak solution to problem 

Proof. Arguing as in Lemma|4l if limfe_i,oo sup^^jQ Wi(A^, /i^) = we know that, up to subse- 
quences, 

t t 
J (jydA'l ^ J (t>dnt, J J v[A'^] ■ V <f) dA'; dr ^ J J v[iir] ■ ^ 4> d^ir dr (fc ^ oo). 

Rti Rti Rd 

The claim of the theorem then follows from Eq. ([^5)) and the arbitrariness of (j), provided we have 
also /jgji^dAo — > J-f^a't'dp- when fc — >■ oo. To show the latter fact, it is enough to prove that 
Wi{Xq, p.) ^ when /c -> oo. 
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Recall that dA§ = dx, with Pq given by Eq. (PO]) . Fixing ip £ Lipi{M.'^) and reasoning like in 
the proof of Lemma [9]- (0) , we find that 

J (f dAg ~ y ^ "^Mj 

where s{x) = X^ieZ'' ail£fc(x) and ai = /gfc ip{x)dx. Therefore: 

y ipd{fi- \l) ^ j {ip- s)dp^^ j {ip{x) - ai) dfl{x). 

But |iy9(a;) — ail < /gt |</3(a;) — <p(j/)| dy < J^k \x — y\ dy, so that from the previous calcula- 
tions we deduce 

J d{fi - X^) < ^ J j \x-y\ dydfl{x) < Vdhk, 

whence Wi{\q, p) < \/dhk — s- for fc oo as desired. □ 

We conclude this section with a simple criterion which implies the convergence property as- 
sumed in Theorem [T^ 

Corollary 13. Assume there exists a bounded set K cM."^ such that supp X'^ C K for all < n < 
Nk and all k > Q. Then (A^)fe>o converges in C([0, T]; Vi{M.'^)) to a weak solution of problem ([¥]) 
when k ^ CO. 

Proof Fix fc > 0, i e [0, T] and let n be such that t^^ < t < i^^^. Let moreover A e B(R'') be 
contained in A'^, then (A) = for all i, k and therefore 

A'AA) = (^1 - ^) X'^iA) + ^A^+i(A) = 0. 

From the arbitrariness of fc, t it follows suppAj^ C K for all fc > and all t € [0, T]. Since K is 
bounded, we can find a ball i?i?,(0) with radius i? > so large that K C 5^(0). Consequently 

\xf dAf(x) = J \xf dA'^ix) <RP < +00, 

i.e., the Aj"s have uniformly bounded moments of any order p > 0. This implies that {A('}fc>o is 
relatively compact in 'Pi{M.'^) for all t E [0, T], which, together with the equicontinuity of the family 
{A^}fe>o entailed by the estimate (|22p . allows us to apply Ascoli-Arzela's Theorem and conclude 
that the sequence (A^')fe>o is relatively compact in C([0, T]; 'Pi{M.'^)). Thus, up to subsequences, 
it converges to some /i, £ C{[0, T]; 7'i(R'^)), which is a weak solution to problem (|4|) because of 
Theorem [H □ 

Remark (CFL condition). Let us define ak ■= V/Pk for all fc > 0. Then the condition expressed 
by Assumption [3] is equivalent to 

VAtk^akhk, (26) 

which implies |7^(a;) — a;| < akhk all x e M'', i.e., the displacements produced by the mapping 
7^ are bounded by the quantity ctkhk- As a consequence, the number of nonempty intersections 
■j^{Ej) n E'f , J S Z'', is finite and bounded from above uniformly in i, which in particular makes 
the series in Eq. (|19p actually a finite sum for all i. 

We observe that Eq. (pS)) is a generalization of the Courant-Friedrichs-Lewy (CFL) condition, 
allowing ak — >■ -l-oo when fc — > oo to ensure convergence to continuous-in-time-and-space solutions 
as stated by Theorem[T21 With the time step At frozen, Eq. ([25)) is used in j5] to prove the stability 
of the numerical scheme (|19p in approximating the solutions to the discrete-in-time model ([8]). 
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5 Models of crowd and swarm dynamics 

In [SJ [5] a class of models describing the collective dynamics of swarms and crowds has been 
introduced. These models are based on the idea that each individual of the group is an intelligent 
(or active) agent, able to develop a behavioral strategy for pursuing specific goals. Since agents 
are not passively dragged by external forces, the Newtonian approach is abandoned in favor of a 
kinematic one, in which the velocity of the agents stems from few basic behavioral rules. These 
concepts are formalized in [S], where it is shown that the above-cited models can be given a 
common formulation within the framework of Eq. (|3]). In particular, the corresponding velocity of 
the agents is 



Adopting the terminology of the kinetic theory for active particles, see [5], we call test agent an 
agent potentially concerned with interactions, whose position is described in the above formula by 
the variable x, and field agents the agents distributed in space (variable y), which the test agent 
may interact with. In more detail: 

(i) Wd : K'' — > K'^ is the test agent's desired velocity, i.e., the velocity of an isolated test agent in 
the absence of interactions; 

(ii) Ux is the interaction neighborhood of the test agent. It conveys the idea that the test agent 
experiences nonlocal interactions with some selected field agents, namely those inside Ux- 

(iii) : K'' — > K is a cut-off function for the set A C M'^, such that 



In particular, the effect of the term xc/x in Eq. (j27|) is to rule out interactions of the test 
agent with field agents outside Ux', 

(iv) r : — > R*^ is the direction of the interaction, depending on the relative position of the 
interacting agents and oriented in such a way that r{y — x) ■ {y x) > 0, with in addition 



(v) / : [0, -l-oo) — >■ K is the interaction strength, which depends on the distance between the 
interacting agents. If / < the interaction is repulsive, i.e., the test agent tries to avoid 
local aggregation with the field agents (as it is common in crowds under normal - i.e., non- 
panic - conditions). If / > the interaction is attractive, i.e., the test agent aims at staying 
close to the surrounding field agents (as in swarms, where group cohesion is advantageous 
e.g., for food search or predator avoidance). 

The shape of Ux depends partly on the criteria used for selecting the field agents to interact 
with. In general, when Ux is symmetric w.r.t. x, interactions are said to be isotropic, as opposed 
to anisotropic when Ux is not symmetric. Notice also that the interaction integral in Eq. ()27[) is 
actually computed w.r.t. the mass measure Nfit, because interactions depend on the number of 
field agents in Ux- 

Prototypes of the above objects are, among others: Ux — Bji{x) for isotropic interactions, 
Ux = S'^{x) for anisotropic interactions (5*^(2;) being the sector of Br{x) with angular width 
a G (0, 27r]), X(7x = 1(7=. > /(|y ~ ^1) = ~1/ \y - A for repulsive interactions, f{\y ~ x\) = \y - x\ 
for attractive interactions, r{y — x) = (y — x)/ \y — x\ (i.e., the unit vector in the direction oiy — x)- 
Although meaningful from the modeling point of view, some of these choices need to be adapted 
in order for the velocity field ^7} to comply with Assumption [TJ 

Concerning this, let us introduce the vector-valued function F : M'* — > R**, 




(27) 



Xa{x)=0 Va;eA^ 



F{x) ■-= f{\x\)r{x). 
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so that (P7)) becomes 

v[fxt]{x) =Vd{x)+N I F{y - x)xu^ (y) dfitiy). (28) 



Such a form encompasses a broader class of models of coordinated behavior based on Eq. ([3|) , for 
instance those in [i]. From now on, we will take (|28p as our velocity model. 

Assumption 4 (Properties of v as in Eq. ([25)) ). We assume that the velocity field in Eq. (|28p 
satisfies the following properties. 

(i) fd is Lipschitz continuous and bounded in W^. 

(ii) There exists R > such that, for each x G Mf^, the set Ux C Bii{x) is measurable and 
isometric to a reference set ?7o ^ -B_r(0). 

(iii) F is Lipschitz continuous in i?^(0). 

(iv) XA is Lipschitz continuous in M'' for each measurable A C R''. 
Remark. Some comments on Assumption |4] are in order. 

(i) In the sequel we denote 

Vd := sup |wd(a;)| < +oo. 

(ii) Saying that Ux is isometric to Uq means that there is a rigid transformation : R'' ^ M'' 
mapping Uq onto Ux'- S,x{Uo) = Ux- Specifically, ^x has the form 

Cx{z)=nxz + x, (29) 

where Tlx G R*^^** is a rotation matrix possibly depending on the point x. 

(iii) Lipschitz continuity of F is required only in the ball i?j^(0), not in the whole space, i.e., 
the condition \F{z2) ~ F{zi)\ < Lip{F) \z2 — zi\ must hold only for \zi\ , \z2\ < R- This 
implies that F is bounded in Bb{0), that is, |-F(2:)| < Cp for \z\ < i?, with Cp ■= 
max{|F(0)|, Lip(F)}(l + i?). 

(iv) The function XA can be thought of as a mollification of 1^. In particular, by continuity 
XA = on dA. Because of the isometry between Uq and Ux, one can check that the 
following relation holds: 

xuAy) = xuoiC\y)). VyeR'^. (30) 

In the next two sections we show that, with both isotropic and anisotropic interactions. As- 
sumption m is sufficient to apply our existence and approximation results to the above-mentioned 
crowd and swarm models. First, we consider the case of a spherical neighborhood, which is the 
most important prototype for isotropic interactions. Later, we extend the analysis to any bounded 
neighborhood, including the significant cases of anisotropic interactions. 

5.1 Spherical interaction neighborhood 

In case of isotropic interactions we set Ux — Bji{x), and consequently Uq = Bii{0). Since Br[x) 
is invariant under rotations, the mapping S^x is simply i,x{z) = z + x. 

By the change of variables z = £,x^{y) in Eq. (pS)) . and recalling furthermore Eq. ([50)). we can 
rewrite the velocity as 

v[^h]{x) = Vd{x) + N J F(z)xs„(o)Wc«(C'#Mt)W, 

where C^#Ht is in Pi{R'^) for all /x* G Pi(R'^) and all fixed x. 

We state a preliminary result, which is useful to verify that such a velocity field complies with 
Assumption [T] 
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Lemma 14. The function z i— > F{z)xBn(o){z) is Lipschitz continuous in M.'^. 



Proof. To study the expression 



e(zi, Z2) := \F{z2)xBR{0)iz2) - F {z-y)xb R{id]{z\)\ , zi, Z2 € M' 



d 



it is convenient to distinguish three cases, 
(i) If zi, Z2 G -B_r(0) we have 

e{zi, Z2) < \F{z2)xBn(o)iz2) - F{z2)xb„Xo){zi)\ + |^(^2)Xi3B(o)(^i) - F{zi)xbr{o){zi)\ 
< (CFLip(xB„(o)) + Lip(-F^)) 1^2 - 21 1 . 

(u) If zi ^ Bfi(0) and Z2 G -Bh(O) (or vice versa) then e{zi, Z2) — \F{z2)xbr(o){z2)\- Let 
zg := Ozi + (1 — 0)z2^ 9 e [0, 1], be a point of the segment connecting zi to Z2 and pick 9 
such that Zfj G 9i?j^(0). Since XBh(o)(%) = Oj have 

e(zi, Z2) = |-F(^2)xsb(o)(^2) - F{zg)xBR(Q){zs)\ < (Cf Lip(xBfl(o)) + Lip(-F')) 1^2 - zg\ . 

On the other hand. \z2 — zg\ = 9 \z2 — zi \ < \z2 — Zi\. 

(iii) If zi, Z2 ^ Sfl(O) then e(zi, Z2) = < |z2 - zi]. □ 

Proposition 15 (Velocity with isotropic interactions). Let Assumption^ hold with Uq — Bji{Q). 
Then the velocity field (|28p complies with Assumption{J[ 

Proof. In the sequel, a; g R'' and /i e 'Pi{M.'^) are fixed but arbitrary. Wc verify the items of 
Assumption [T] in order. 

(i) The field v is uniformly bounded thanks to 



\v[ti]ix)\ <Vd + N \Fiz)xBRio)iz)\ ci(4-i#/i)(z) <Vd + NCpl^iU^")) = + NCp- 




(ii) As for Lipschitz continuity, we check it separately w.r.t. to x and to /x. 
(ii-a) Let xi, X2 E M'', then, recalling Lemma UM 



\v[fi]{x2) - v[^i]{xi)\ < \vdix2) ~ Vdixi)\ + N / F(z)xs«(o)(^)d(C#M-C#M)W 




< Lip(«d) \X2 - xi\ + NUpiFxBR{o))Wii^-^^#fi, C/#Ai). 



In addition, 




sup 

ipeLipi(R<i) 





whence the Lipschitz continuity of a; i— >■ v[^]{x). 
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(ii-b) Let now jJ., v £ ViiW^), then, invoking again Lemma [TH 



\v[u]{x)-v[^J]{x)\^N 



F{z)xB,,(0){z)d{£,-^#j^ - C'#m)(^) 



< iVLip(Fxi3«(o))VKi(C'#M, 
= iVLip(Fxi3„(o)) sup I LP o d{v - ^i) 



'/=eLipi( 

^ NU^(FxBn[0))Wl{^i, v), 

since € Lip]^(R'^) implies that (po^^^ spans the whole space Lip2(M'*) when (p varies 
in Lip]^(K'^). Hence also [l i-> is Lipschitz continuous. 

(iii) Finally we examine the linearity w.r.t. to the probability for convex combinations. Let 

II, V £ 7'i(M'') and a e [0, 1], then 



v{a[i ^ {Jl - a)v\[x) ^ va{x) + aN J F{z)xBn(o)iz) d{C^^#^){z) 
+ {l-a)N f F{z)xB^^o){z)d{C'#i^){z)- 



Writing Vd{x) = avd{x) + (1 — a)vdix) and collecting the coefficients a and 1 — a gives the 
result. □ 

Owing to Proposition we can state that models based on Eq. ^ , with a velocity field 
featuring isotropic interactions as in Eq. (j28|) . see e.g., [8], admit probability measure solutions for 
any initial distribution of the agents with finite first and second order moments (e.g., an initial 
distribution with compact support). Moreover, such solutions can be duly approximated using 
the numerical scheme ([T^. 

Remark (Unbounded interaction neighborhood). As a modification to Assumption |4l we can allow 
Uq = Ux = M'' along with the boundedness assumption F{x) < Cp for all x £ R''. The above 
arguments can be promptly adapted to show that also in this case the velocity field ([^5)) complies 
with Assumption [TJ However, it should be noted that, for the applications we have in mind, the 
physically relevant cases are those in which the interaction neighborhood is bounded. 

5.2 Bounded interaction neighborhood 

In this section we drop the specific hypothesis that Ux be a ball. We allow it to have a generic 
shape, with the only requirement of being bounded, i.e., contained in a ball. This encompasses 
the important case of anisotropic interactions. In such cases, the neighborhood Ux may not be 
invariant under rotations, therefore we have to consider the full form (j29|) of the transformation 
^x- Performing again the change of variables z = £,x^{y) in the integral the velocity takes 

now the form 

vifitKx) = v^x) + N f F{nxz)xuo{z)d{£,-^4h^it){z) 



with the rotation matrix explicitly appearing in the argument of the function F. To deal with 
it, for the sake of simplicity we confine ourselves to the two-dimensional case (d = 2), for then a 
simple representation of TZx is available: 



Tlx 



cos "dx — sin dx 
sin "dx cos "dx 
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•dx G [0, 2tt) being the angle of rotation which determines the local orientation in the point x of 
the neighborhood of interaction. 

The choice of "dx has to do with the way in which the anisotropy of the interactions is modeled. 
In the models of crowd and swarm dynamics we are considering, is the angle formed by the 
desired velocity V(i{x) w.r.t. a fixed reference direction, for instance the horizontal one identified 
by the unit vector i. Assuming for simplicity that Vd has constant unit modulus, this implies 



In the second formula, Vd{x) and i arc thought of as embedded into M.^, with x denoting vector 
product and k the unit vector orthogonal to the plane of V(i{x) and i. 

Remark. More in general, Eq. ([3T|) holds with Vd{x) replaced by Vd{x)/ \ vd{x)\, which is a Lipschitz 
continuous field provided |wd| is uniformly bounded away from zero. 

The formalism introduced above allows us to prove the following technical fact. 
Lemma 16. For all xi, X2, z G it results 



The same holds with TZxi, 'R-X2 replaced by Ti-xl , Ti-x} ■ 
Proof. A straightforward computation shows that 

\{Tlx2 " 'R-xi)z\ = V - cos-dx^f + (sin??^;^ - sindxiY \z\ , 

and the same is true also using inverse matrices. In addition, 

\cOS'dx2- COST^xA = \{vd{x2)-Vd{xi))-\\<\vd{x2)-Vd{xi)\, 

\sv[i'dx2-sv[i'dxx\ = |[(ud(a;2) - Wd(a;i)) X i] • k| < |wd(a;2) - Wd(a;i)| , 

hence the thesis follows from the Lipschitz continuity of v^. □ 

With Lemma[Tn]we are in a position to prove that the velocity ([25)1 complies with Assumption[T] 
also in case of anisotropic interactions. 

Proposition 17 (Velocity with anisotropic interactions). Let Assumption^ hold and assume 
d = 2. Then the velocity field (|28p complies with Assumptions^ 

Proof. (i) Uniform boundedness, Lipschitz continuity w.r.t. the probability, and linearity w.r.t. 
to the probability for convex combinations follow straightforwardly from calculations entirely 
analogous to those performed in the proof of Proposition llSI In fact, it is sufficient to observe 
that, TZx being an isometry, the function F{TZx-) is Lipschitz continuous and bounded in 
Br{0) with Lip(i^(7?.2--)) = Lip(^) for all x £ W^. Moreover, by the same argument as in 
the proof of Lemma [HI the function F{Tlx-)xUo{') is Lipschitz continuous in M'* with the 
same Lipschitz constant as Fxuoi thus in particular independent of x. 

(ii) Lipschitz continuity w.r.t. to x is instead more delicate, because it involves directly the 
rotation matrix TZx. Let si, X2 € and fix € 7'i(R^), then 



cos-dx = Vd(x) ■ i. 



sini?^ = {vd{x) X i) • k. 



(31) 



\{nx2-nx,)z\ < V2Uj>{vd)\x2~xi\\z\. 



\v[fi]{x2) - v[fi]{xi)\ < \vdix2) - Vd{xi)\ 




<Up{vd)\x2~xi\ + N / |F(7^,,z)-F(7^,,z)|xc/o(^)c^(C/#M)(^) 




(32) 
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In the first integral at the right-hand side of ([5^ we can assume z E Uq, for otherwise 
XUoi^) = 0- Hence \z\ < R and moreover TZxZ G -Bfl(O) for all x G . Lipschitz continuity 
of F in that ball, along with Lemma [TBI implies 

|F(7^,,2)-F(7^,,z)| < U^{F)\{Tlx,-Tlx,)z\ < V2 Lip(F) Lip(wd)i? | X2 — Xi \ , 

so that finally 

N [ |F(7^,,z)-F(7^,,z)|xc/oWd(C.'#/^)W<^^ / |F(7^,,z) - F(7^,,z)| d(C,¥A*)W 



< NV2Up{F) Up{vd)R \x2 - xi I , (33) 
having observed that (CK2^#M)(t^o) = ^^{Ux2) < 1- 

As far as the second integral at the right-hand side of ([5^ is concerned, we have 
N 



F{nx,z)xuo{z) - C/#m)(^) 

<N j |F(7^,,c.'(y))x[/o(C(^/))-^^(^.lC'(y))xao(C'(y))| ^^(y) 
= iv / |i^(7^.,C/(y))xao(C/(y))-^(7^.,C(2/))xc/o(C(2/))|rfM(2/)• (34) 



Notice that we can confine ourselves to y E UxiUUx2 , for otherwise xua iCx ^ (j/)) ^ XU:^ . (y) = 
for both j = 1, 2. We distinguish two cases. 

(ii-a) \x2 -xi\> 2R. 

In this case Ux^ n Ux2 = because the balls Bji{xi), Br[x2) are disjoint. Thus: 

m^N I \F{TZx,cHy))\xu.Ay)dfxiy) + N J |F(7^,,c/(y))| Xa., (2/) ^^(y)- 

For aU y G C/,^., j = 1, 2, it results ^"/(y) G C/q C ^^(0), hence 7^,,C/(2;) G Br{0) 
and we can use the boundedness of F in that ball to get 



< NCf 



J xu., (y) dy.{y) + j xv.., (y) d^i{y) \ < NCp^U.^ U C/, J < NCp. 



But 1 < 2R ' therefore we conclude 

m<^\^2-x,\. (35) 

(ii-b) - xil < 2R. 

In this case the neighborhoods Ux^ , C/a;2 need not be disjoint, however we can resort to 
the Lipschitz continuity of the function F{TZx-)xUo{')'- 

m < Lip(^Xc/o )^ / (y) - C (y) I dti{y) 

- Lip(Fxc/o)^ J \n-^\y - X2) - n-^iy - xi)\ dfi{y) 



Uv{FxUo)N j \{nxl - UxDiy - xi) + 7^-^l(x2 - xi)\ dn{y) 

t/x, U(7x, 
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and further, thanks to Lemma [T51 and to the fact that TZ^^ is an isometry. 

/ \ 

<Up{Fxuo)N V2Up{vd) J \y - xi\ dn{y) + fi{U,, U U,,) 



\X2 - Xi\ 



Let us examine the term with the integral. If y G Ux-^ then |?/ — a;i| < R whereas 
if y e Ux^ then — Xi| < \y — X2 \ + \x2 — xi\ < 3R. Finally, |j/ — xi| < 3R for all 
y eUx^U Ux2, which says 

dMl) <Lip(^^Xc/o)^(3\/2Lip(t;d)i?+l)|x2-xi|. (36) 



In conclusion, from ([35)) and (j36|) we deduce that there exists a constant C > such that 
< C\x2 — xi \ for all xi, X2 G K^- This, together with the estimate ([55)1 . completes the 
proof of Lipschitz continuity of the mapping x i— >■ v[^]{x). □ 

In view of Proposition [T7] we conclude that two-dimensional models based on the velocity 
(|28p with anisotropic interactions have probability measure solutions, which can be approximated 
arbitrarily well using the scheme (jl9[) on finer and finer numerical grids. Notice that, as far as 
e.g., crowd dynamics is concerned, two-dimensional problems are enough for applications. 

Remark (Higher dimension). For d > 2 additional technicalities arise, due to a more complex 
structure of the rotation matrix. Nevertheless, in the special case that the desired velocity is 
constant in x, it is straightforward to extend the results to any spatial dimension. In fact, the 
rotation matrix being independent of x, Proposition [T7] can be proved without using Lemma |16l 
which is the only point where we use the explicit representation of the matrix. Models with 
constant desired velocity have been recently proposed for swarm dynamics problems [5] and for 
rendez-vous algorithms [3]. 

Remark (Zero desired velocity). When interactions are anisotropic and the desired velocity is zero 
[5], the orientation of the neighborhood of interaction cannot be defined by Eq. ([5T|) . However, 
this issue can be solved by replacing Vd in pip with any other Lipschitz continuous unit vector 
field, e.g., a nonzero constant one, with the only purpose of defining a rotation angle. Clearly, this 
problem does not arise if the desired velocity is zero but interactions are isotropic, as in 



6 Case study: discrete models 

In this last section we put the theory into operation by giving an example of explicit solution 
to the Cauchy problem Furthermore, wc visualize the convergence to such solution of the 
approximations produced by the numerical scheme discussed in Section |4l 

For illustrative purposes, we consider the simple case of a system of agents featuring isotropic 
interactions with zero desired velocity: 

v[fit]ix) = N J F{y- x)xBR{x){y)d.Ht{y). (37) 

Moreover, we prescribe as initial condition the discrete probability measure 

1 ^ 

1=1 

where (5^i^ is the Dirac mass concentrated in Xq (i.e., for any A € B{M.^) it results S^i^{A) = 1 if 
Xq G a, Sxi.^{A) = otherwise) and x^, . . . , Xq are N selected points in M'*. 

We recall that jl is the common law of the random variables Xq, i = I, . . . , N, expressing 
the initial positions of the agents. The structure of p, implies that each Xq is a random 



28 



variable taking almost sm'ely the values xj, , each with probability 1/N. Therefore we 

are considering a discrete model of a group of indistinguishable agents initially concentrated in 
Xg, . . . , Xq . The indistinguishability is reflected by the fact that any of the points Xq can be, with 
equal probability, the initial position of the generic i-th agent. 

We find a solution to the Cauchy problem with initial condition ([551) by the method of the 
characteristics (cf. Section . In particular, we know from Eq. ^ that /if — 7t#M; where jt is 
the flow map. From the linearity of the push forward we first deduce nt = lt4l^^x' j then 



otherwise. 



we observe that for any measurable set A it results 

Hence jt^S^i^ = S^^^^i^-^ and we can write the solution as 

^E'^^.K)- (39) 



N 

i 



N 
1=1 

Recalling that /ij is the law of the random variables X'j: ~ K[XI\Xq], i = 1, . . . , N, from Eq. (j39|) 
we infer that each XI takes almost surely the values 7t(xJ), . . . , Jtixo), each with probability 1/iV. 
Therefore, at every time t > the distribution of the group of agents is concentrated on the discrete 
set of points 7t(xJ), . . . , 7t(x^). Notice that again we cannot associate deterministically a given 
agent with its position because of the indistinguishability of the agents. However, we can describe 
the trajectories of the agents by means of the mappings 1 1-> 74 (xg), 1 = 1, ■ ■ ■ , N. 

The flow map is defined by Eq. ([7]), which with the velocity (|37)) and the initial condition ((38|) 
yields 

= E^(^*(4) " ^t(^))XBuht(x))Mxi)) ^^^^ 
^70(2;) = X 

for all X € M'^. The discrete structure (p9| of /.it makes it actually sufficient to solve problem (|40 



for X = Xq. Setting xi{t) := 7t(a;Q) and computing Eq. (|40p for x = Xq we find the initial- value 



problem 

N 

-.-E^(-.--OXB«(.)(-.) (/ = i,...,^), (41) 

xi{0) = x[, 

thus we conclude that constructing the measure ([5^ amounts to solving the system of ODEs (HU, 
whose solutions are the positions of the agents at t > 0. 

From the numerical point of view, we can either approximate the measure (|39|) by using the 
scheme p9p or integrate problem l|4ip via a standard numerical method for ODEs. The remaining 
part of this section is devoted to show the convergence of the approximations obtained from (jl9p 
to the numerical solutions of (j4ip in a toy model. Let us consider a one-dimensional problem 
(d = 1) with = 10 agents, whose initial positions are sampled from a uniform distribution on 
the interval [0, 1]. Agents repel each other according to the following repulsion function: 



max^{N,e} ^"-^ > 



which is from the product of 



max{|z|,e}' max{|z| , e} 

The repulsion strength / is inversely proportional to the distance between the interacting agents 
(up to a minimal threshold e) and the direction of the interaction r is a Lipschitz mollification of 
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t = T/2 



t = T 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Figure 1: Numerical solution to Eq. ([3]) with discrete initial datum ([55]) at two successive times. 
Bullets are the solution of the ODE system (PT|) computed by the explicit Euler scheme. The 
continuous curve is the solution computed by the numerical scheme (|19p on meshes with different 
level of refinement. 



the unit vector z/ \z\. The reference interaction neighborhood is Bjf {0) — {—R, R) with cut-off 
function ^ 

XBH(0)(^) = e~^?^lB,(0)(z) (6>0), 

which is a C°° mollification of the indicator function of Bii{0). Parameters are set to R ~ 0.1, 
a = 0.01, b = 0.02, e = R/A. The computational time is T 0.1. 

Simulations of the ODE system (pij) and of problem (U) with initial condition were run 
independently and their results visualized on the same graphs of Fig. [1] for duly comparison. In 
particular, the ODE system was numerically integrated using an explicit-in-time Euler scheme, 
whereas the conservation law for the probability fit was solved through the scheme (fT9)) on different 
meshes, choosing 




(0 < (5 < 1). 



Note that this entails l3k ^ h]^^ in Assumption [3] and ak ^ h^i^^ hi Eq. (pS)) . Figure [T] displays the 
numerical solution computed with 5 = 0.9 and /i^ = 1/A: in the three cases k = 10^, 10"^, 10**, at 
two different time instants. Convergence toward the exact solution p9|) as the mesh refines can be 
visually appreciated, although approximating singular measures with densities is a difficult task, 
which requires very fine and computationally expensive meshes to get accurate results. Therefore, 
when such a structure of the solution is numerically sought, it is more efficient to exploit the stated 
equivalence of the original problem with the discrete system of ODEs. 
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